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Abstract 

Diffractive parton distribution functions give the probability to find a parton 
in a hadron if the hadron is diffractively scattered. We provide an operator 
definition of these functions and discuss their relation to diffractive deeply 
inelastic scattering and to photoproduction of jets at HERA. We perform 
a calculation in the style of "constituent counting rules" for the behavior of 
these functions when the detected parton carries almost all of the longitudinal 
momentum transferred from the scattered hadron. 
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I. INTRODUCTION 



Recently, the Zeus and HI experiments at HERA have reported the first evidence for 
diffractive deeply inelastic electron scattering [|TJ, 

e + A^e + A' + X. (1) 

This is an example of a more general phenomenon, diffractive hard scattering, in which a high 
energy incident hadron participates in a hard interaction, involving very large momentum 
transfers, but nevertheless the hadron itself is diffractively scattered, emerging with a small 
transverse momentum and the loss of a rather small fraction of its longitudinal momentum. 
One may say that the hadron has exchanged a pomeron with the rest of the particles 
involved and that the pomeron has participated in the hard interaction. The possibility 
of such interactions was proposed by Ingelman and Schlein || on the grounds that the 
entity exchanged in elastic scattering, called the pomeron, must be made of quarks and 
gluons, which, being pointlike, can participate in hard interactions. The theoretical ideas 
and formulas involved are elaborated in some detail in Ref. ||. The predicted phenomenon 
was seen in jet production in hadron collisions by the UA8 Collaboration Q. 

As discussed in our previous work ]5[], the Ingelman-Schlein model can be thought of 
as involving a "diffractive parton distribution function," which is the subject of this paper. 
The idea is that this function, 

dxjp dt 

represents, in a hadron of type A, the probability per unit <i£ to find a parton of type a 
carrying momentum fraction £, while leaving hadron A intact except for the momentum 
transfer characterized by parameters [xjp,t). Here t is the invariant momentum transfer 
t = (P A - P A ,) 2 while xjp is the fraction of its original longitudinal momentum lost by the 
hadron. The parameter /j is the factorization scale, roughly, the resolution of the parton 
probe. A function expressing the same physics as the diffractive parton distribution (0) has 
been proposed by Veneziano and Trentadue under the name of "fracture function." The 
details are a little different, as we will explain in Sec. JTI]. The original paper of Ingelman 
and Schlein did not mention the function (0) but instead introduced a related function, the 
"distribution of partons in the pomeron." 

Our purpose in this paper is, first of all, to relate these various functions and the ideas 
behind them to one another and to comment on the likely validity of the formulas that 
express cross sections in terms of these functions. We give operator definitions for the 
diffractive parton distribution functions and discuss the evolution equation that they obey. 
We briefly review the expected behavior of /x; xjp,t)/dxjpdt for small /?, where (3 = 

£/x]p. Then we use a perturbative calculation to explore how these functions behave for 
small values of 1 — (3. In the Ingelman-Schlein language, our results favor a rather "hard" 
distribution of partons in the pomeron. We also present, in an appendix, a calculation of 
d fa/A(^ fJ>',Xjp,t)/dxjpdt in a simple model. We conclude with some observations on the 
experimental consequences of the theory. 
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II. DIFFRACTIVE DEEPLY INELASTIC SCATTERING 



In a deeply inelastic scattering reaction, a hadron A with momentum Pj{ is struck by 
a far off shell photon with momentum q^. It is convenient to use momentum components 
k^ = (k + , k~ , k), where fc ± = 2 -1 / 2 (/c° ± k 3 ), and where we denote transverse components 
of vectors by boldface. We work in the brick wall frame, in which P% = (P^, M^/[2P^], 0) 
and = 2~ 1 / 2 (—Q, Q, 0). One measures the standard hard scattering variables Q 2 = —q ■ q 
and x = Q 2 /[2Pa ■ q\. In some deeply inelastic scattering events there will be in the final 
state a diffractively scattered hadron A' with momentum 



P A , = [1 - x p \P2 



n + mi 

2[1 - xjp]P^ 



+ ' 



(3) 



as in Fig. (|]). The hadron has lost a fraction xjp of its plus momentum and has gained 
transverse momentum P^/. The invariant momentum transfer from the proton, t = (Pa — 
Pa>)\ is 



p2 

" A' 



+ x 2 P M 2 



I -Xjp 



(4) 



The events in which we are interested have small t. One expects \t\ ^ 1 GeV to be 
typical. We also suppose that xjp is rather small. One expects pomeron physics to be 
dominant for xjp < 0.1. Having found such events, one can construct the contribution to 
F2 from final states containing a diffractively scattered hadron with variables t and xjp: 
dFf s (x, Q 2 ; xjp, t) jdxjp dt. 




FIG. 1. A graph for e + p — > p + X. 



The model of Ingelman and Schlein as applied to deeply inelastic scattering is simple 
to state. We begin with the usual factorization theorem for the structure function F2. 



F 2 (x,Q 2 )=J2 J\ f a/ A(Z,li) F 2 ^x/^Q 2 -fi). 



(5) 



Here f a /A{£,, A*) is the distribution of partons of kind a in hadron A as a function of momentum 
fraction £, as measured at a factorization scale /1, while F 2 a is the structure function for 
deeply inelastic scattering on parton a. If, for simplicity, we ignore Z exchange, then F 2 is 
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F 2 , a (x/£, Q 2 ; fi) = e 2 a 5(1 - x/£) + 0(a s ). 



(6) 



Thus F 2 is rather trivially related to the parton distribution functions at the Born level; 
nevertheless, conceptually the distinction between F 2 (x, Q 2 ) and fa/A^,^) is quite impor- 
tant. As in our previous paper ||, we break the analysis into two stages. In the first stage, 
we hypothesize that the diffractive structure function F^ 15 can be written in terms of a 
diffractive parton distribution, Eq. (0): 

2^ L j„ m F 2A {x/Z,Q (7) 



dxjp dt ^ Jo dxjp dt 

In the second stage, we hypothesize that d f^(x a , fi) / dxjp dt has a particular form: 

dxjp dt = 8^2 x ip fa/AU^IP, t, fi) . (8) 

Here (3a (t) is the pomeron coupling to hadron A and a(t) is the pomeron trajectory. We 
distinguish the "Regge factorization" ofEq. (j^) from the "diffractive factorization" ofEq. (0). 

In Eq. we adopt standard conventions such that the proton-proton elastic scattering 
amplitude is 

M = -(3 p {t) 2 s a(t) . (9) 
Then the elastic scattering cross section is 

| = J_ |/3p(i)|V M'M». (10) 

while the total proton-proton cross section is 

ff«(BP) = ^(O) 2 ] . (11) 

The normalization factor l/(87r 2 ) in Eq. (|8|) is quite arbitrary. Here, we have adopted the 
convention of Donnachie and Landshoff 0. 

The function f a /ip(P, t, fi) thus defined is the "distribution of partons in the pomeron." In 
writing Eq. (|8|), one thinks of the pomeron as a continuation in the angular momentum plane 
of a set of hadron states. Since hadrons contain partons, the pomeron should also. Thus 
one has in Eq. (§) the standard factors describing the coupling of the pomeron to hadron A, 
together with a distribution of partons in the pomeron |2]||. Inserting Eq. (|[) into (0), one 
obtains the model of Ingelman and Schlein, applied to the case of deeply inelastic scattering: 

dF 2 diff (x,Q 2 ;W) \P A (t)\ 2 i-M^f 1 



dxjp dt 8n 2 



%T 2Q(t) W df3f a/1P 0,t,fi)F 2 , a (f3/(3,Q 2 -,fi), (12) 



where (3 = x/xjp. We offer here a word of caution. Both the structure of Eq. fll2"D and 
the language "distribution of partons in the pomeron" suggest that the hadron emits a 
pomeron some long time before the hard interaction and that the pomeron then splits into 
partons, one of which participates in the hard interaction. This interpretation is, however, 
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not required by Eq. ( pT2| ) and is surely quite misleading. In a diagrammatic interpretation 
of pomeron exchange ||, the exchanged quanta have small plus and minus components of 
momentum. Thus the exchange takes place over a long interval (Ax + , Ax~) in space-time. 
It begins long before the hard interaction and ends long afterwards. Our diagrammatic 
analysis in Sees. |V1I| and |V111| will provide an illustration of this picture. 



We see that Eq. (0), can be regarded as a version of the Ingelman-Schlein model (|12|) that 
is more parsimonious in its assumptions. Eq. (^) says only that factorization still applies 
when hadron A is diffractively scattered. The Ingelman-Schlein model flT2|) assumes that 
Regge phenomenology is applicable and, with the aid of this assumption, has more predictive 
power. 

In this paper, we concentrate on the case in which hadron A' is the same kind of hadron 
as hadron A, so that vacuum quantum numbers are exchanged, and we consider xjp to be 
small enough so that pomeron exchange dominates. One should keep in mind, however, that 
Eq. (0) admits generalizations to cases where A' ^ A and where xjp is not at all small. One 
can also generalize Eq. (|T2|) to A 1 ^ A, but then xjp should be fairly small in order that just 
one or two Regge exchanges dominate. 

The diffractive factorization equation (|7]), or rather a very closely related equation, 
has been introduced by Veneziano and Trentadue ||. These authors call the analogue 
of d fa/A(^ fj) / dxjp dt a "fracture function." Stated precisely, a fracture function is 

1ZT-T+ • ( 13 ) 



dxjp Jo dxjp dt 

By integrating over t, Veneziano and Trentadue eliminate a variable that is perhaps of 
secondary importance. However, there is some advantage to not integrating over t. We are 
interested in the physics of diffraction, which occurs in the small t region. If we integrate 
over t, then we are forced to consider also the large t region, in which the hadron A' is to be 
thought of not as the original hadron appearing in a scattered form but as a random hadron 
in a high Pt jet produced in the hard interaction. Taking this possibility into account leads 
to certain complications in the formulas. 

In the following sections, we analyze the diffractive parton distributions. According 
to Eq. (0), the measured quantity dF^/dxjpdt is approximately the sum of diffractive 
quark distributions weighted by the square of the quark charges. There are higher order 
corrections to this relation, some involving the diffractive gluon distribution. Thus these 
distributions are rather directly related to experiment. The reader may wonder why we 
concentrate on the theoretical diffractive parton distributions rather than on the physical 
quantity dF^ jdxjp dt. The reasons are the same as in ordinary hard scattering: 1) the 
diffractive parton distributions are process independent and 2) the factorization ([?]) allows 
one to include perturbative corrections to the hard scattering. The reader may also wonder 
why we don't frame the analysis in terms of the distribution of partons in the pomeron. 
Our excuse is ignorance. We don't know how to relate the Regge factorization in Eq. to 
quantum field theory. 
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III. THE DIFFRACTIVE PARTON DISTRIBUTION 



In this section we give an operator definition of the diffractive parton distribution. We 
write the ordinary distribution of a quark of type j G {u, u, d, d, . . .} in a hadron of type A 
in terms of field operators ip(y + , y~, y) evaluated at y + = 0, y = PJIOl: 



1 1 



si 



f j/A ({, y) = — - £ / dy-e-* p * y " (P A , s A \^(0, y~, 0) 7 + ^(0) \P A ,s A ). 



Similarly, the ordinary distribution of a gluon in a proton is written as 

/ S M(£, A*) = ^pr\ E / dy-e^r {P A , s A \F^(0, y~, 0)^ F a (0) v + \P A , s A ). 



(14) 



(15) 



The proton state \P A , s A ) has spin s A and momentum P A = (P A , M A /[2P A ], 0). We average 
over the spin. Our states are normalized to 



(k\p) = (2tt) 3 2p° 5 3 {p- k) = (2tt) 3 2p + 5(p + - k + ) 5 2 {p - k) 



(16) 



The field ipj(0, y , 0) is the quark field operator modified by multiplication by an exponential 
of a line integral of the vector potential: 



4>j(o,y~,o) 



Vexp(igJ dx~ A+(0,x~,0)t t 



Likewise F a (0, y ,0) +lJ is defined by 



Vexp^ig J dx A^(0,x , 0) t^j 



ipj{0,y-,0) 



F b (0,y-,0) 



[IV 



(17) 



(18) 



The V denotes path ordering of the exponential. The matrices t c in Eq. ([L7]) are the gen- 



erators of the 3 representation of SU(3), while in Eq. (|18 ) they are the generators of the 8 
representation. These operator products are ultraviolet divergent, and are renormalized at 
scale /x using the MS prescription, as described in ||. 

The motivation for these definitions is that in QCD canonically quantized on null surfaces 
x + = const, using A + = gauge, the operators measure the probability to find a quark and 
gluon respectively carrying plus component of momentum equal to £P A . The line integrals of 
the color potential restore gauge invariance. Then MS renormalization removes divergences. 

The line integrals of the color potential have a physical interpretation. Whenever a 
parton is measured by a short distance probe, the color carried by that parton has to go 
somewhere. For instance, in deeply inelastic scattering, the color is carried away by the 
recoiling struck quark. In the definition of the parton distribution function, the recoil color 
flow is idealized as an infinitely narrow jet moving with the speed of light along the path 
x^ = (0, x~ , 0) with y~ < x~ < oo. Any gluons from the color field of the hadron can couple 
to this idealized color source. 

Consider now the diffractive distribution of a quark in a proton. The operator is the same 
as in Eq. fli~4|), but the proton is required to appear in the final state carrying momentum 

P'a- 
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d 3 P A , 

1 1 



4vr2 t , 



J dye-*^- Yl (Pa, s a \^(0, y~ , 0)\P A > , s A >; X) 

X,s A , 

xy + (P AI ,s A r,X\^{0)\P A ,8 A ). (19) 



We sum over the spin s A > of the final state proton and over the states X of any other particles 
that may accompany it. Similarly, the diffractive distribution of gluons in a hadron is 



(2nf2E A , dft f^ ] = G*f A {PA,PA.,t,ri 
d?P A , 

1 1 



2tt£P a + 2 



]T / dy-e-^- £ (P A ,s A \F a (0,y-,O) + »\P A ,,s A r,X) 

sa X,s A , 

x{P A ,,s A ,;X\F a {0)+\P A ,s A ). (20) 



The Green function G„/a f° r a parton of type a can, at least in principle, be computed from 
Feynman diagrams together with the Bethe-Salpeter wave functions for the bound states. 

We will want to change variables to xjp and t as defined in the previous section. Using 
Eqs. (H) and (||), we obtain 

d 3 P A > dxjp 2 „ 1 



2E A 2(1 -xip) 4 
Integrating over the azimuthal angle 0, we have 



d 2 P A , = - dxjp dt d(j) . (21) 



df*fA£,u) i 

^Bt-t^^.p^) (22, 

where G^K for quarks is given in Eq. flTUp and G^ A for gluons is given in Eq. fl2H|). 



IV. EVOLUTION EQUATION 

As mentioned in the previous section, the diffractive parton distributions are ultraviolet 
divergent and require renormalization. It is convenient to perform the renormalization using 
the MS prescription, as discussed in PJIQf. This introduces a renormalization scale [i into 



the functions. In applications, one sets /i to be the same order of magnitude as the hard 
scale of the physical process. 

The renormalization involves ultraviolet divergent subgraphs, such as that shown in 
Fig. 0(a). Subgraphs with more than two external parton legs carrying physical polarization, 
such as that shown in Fig. 0(b), do not have an overall divergence. Thus the divergent 
subgraphs are the same as for the ordinary parton distributions. We conclude that the 
renormalization group equation for the diffractive parton distributions is 

^ dx^dt = ? h T Pa/b{i,Z ' asifi)) dx^dt (23) 
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(a) 




FIG. 2. Renormalization of the diffractive quark distribution. The subgraph in (a) denoted 
by heavy lines is ultraviolet divergent, and thus contributes to the evolution kernel. The subgraph 
in (b) is not ultraviolet divergent (with transverse polarizations for the incoming gluons). The 
diagrammatic notation is that of Ref. || 



with the same DGLAP kernel fill , Pa/b(£,/z, a s (fi)), as one uses for the evolution of ordinary 
parton distribution functions. 

If, following Veneziano and Trentadue, we integrate over t, then the large t integration 
region introduces new ultraviolet divergences and the renormalization group equation is 
modified ||. In this paper, we choose to restrict integrations over t to the small t region. 



V. VALIDITY OF DIFFRACTIVE FACTORIZATION 

In Sec. H|, we presented the hypothesis of diffractive factorization for diffractive deeply 
inelastic scattering, as represented by Eq. (^): 

E /„ <% ^ M F 2)a (x/Z,Q ;//). (24) 



dxjp dt ^ Jo dxjp dt 

This is an example of the more general hypothesis of factorization for other kinds of diffrac- 
tive hard scattering. Another example is diffractive jet production. Consider, for example, 
the inclusive cross section for the production of two jets in a high energy collisions of two 
hadrons, A and B. (At HERA, this would be p + 7 — > jets where hadron B is the the 
hadronic or "resolved" part of the photon.) Let the initial hadron A have momentum 



2 



n = (Pa, PI, Pa) = 7^, 0) , (25) 
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while hadron B enters the scattering with momentum 



M 2 

PS = (Pb, Pb, Pb) = Pg, 0) • (26) 

We specify the two jets by variables Ej<, X A , and Xb, given in terms of the four momenta 
Pf and P<2 of jets 1 and 2 by 

E T = (|P 1 | + |P 2 |), 
X A = (P+ + P+)/P+ , 

X B = (P 1 - + P 2 -)/P B . (27) 

If we add the requirement that hadron A emerge scattered with scattering parameters (xjp, t), 
then we have diffractive jet production. The corresponding hypothesis of diffractive factor- 
ization for this cross section is 

da m (A + B -> A + jets + X) 

■ rsj 

dEx dX A dXs dxjp dt 
^ /' dfffi(x a ,wx P ,t) r da{a + b -> jets + X) 

g I dXa d^dt J dXb fb/B{Xb > A dE T dX A dX B ■ (28) 

Of course, Eqs. fl24|) and fl28|) are approximations, as indicated by the ~ signs. We un- 
derstand the hypothesis of diffractive factorization to mean that the corrections to these 
relations are suppressed by a power of m/Q or tu/Et, where m represents the momentum 
scale of soft hadronic interactions and Q or Et is the scale of the hard interaction. 

This general diffractive factorization hypothesis was put forward by Veneziano and 
Trentadue in their paper on fracture functions. The Ingelman-Schlein model @ de- 
mands diffractive factorization plus the Regge structure of the diffractive parton distribu- 
tions. Thus, in a strict interpretation, the validity of the Ingelman-Schlein model logically 
implies the validity of diffractive factorization. On the other hand, one might interpret the 
Ingelman-Schlein model as being valid if corrections to it, while not vanishing in the limit 
of large Q, were nevertheless numerically small. Thus, for instance, the authors of Ref. |§ 
speculated that the factorization inherent in the Ingelman-Schlein model was not likely to 
be exact up to m/Q corrections but might have the same status as Regge factorization, 
which has proven to be a useful approximation even if it is not exact. 

The hypothesis of diffractive factorization appears to us to be correct in the case of deeply 
inelastic scattering. A detailed proof of this statement is beyond the scope of this paper. 



However, we can briefly sketch how such a proof would go, following the ideas of Refs. fll2| , |13 
The singularities of the cut Feynman graphs for diffractive deeply inelastic scattering are 
such that the leading integration regions involve 1) a beam jet in the direction of the initial 
hadron A (which includes the final state diffracted hadron A'), 2) a hard interaction, 3) one 
or more final state jets that are not in the direction of hadron A, and 4) possible soft gluons 
(sometimes with soft quark loops) that may communicate between the beam jet and the final 
state jets. Factorization would be more or less kinematic were it not for the possibility of 
the soft gluons linking the beam jet with the final state jets. One must use gauge invariance 
to show that the soft gluons don't really "see" the details of the final state jets, so that the 
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connections to the final state jets can be replaced by connections to the idealized jet that 
is embodied in the line integral of the field operator in the definitions of the diffractive 
parton distribution functions, Eqs. ( |14| , p~5| , [T7] , p~8D . 

What of diffractive factorization for processes with two hadrons in the initial state, such 
as p + p —>■ jets or 7 + p — > jets! Here the proof of ordinary factorization is much more 
delicate. The problem is that low momentum gluons can communicate between the partons 
of the two beam jets. This can happen even before the hard scattering takes place, as in 
Fig. |3|. When one looks for the hard process inclusively, such effects cancel JT^JT^]. But 
the cancellation requires a sum over final states. Collins, Frankfurt, and Strikman |TJJ 
have argued that the demand that the final state include a diffractively scattered hadron 
destroys the factorization. In Ref. ||, we looked at this problem in the context of a simple 
perturbative model. We found that the diffractive factorization hypothesis is, indeed, not 
valid. New terms proportional to 5(1 — X^jxjp) appear in Eq. (^). (Presumably in a more 
general model one will also have factorization violating terms that are not proportional to 
5(1 — Xa/xjp) but are singular as (1 — Xa/xjp) — > 0.) These new terms have an interesting 
structure that can be investigated experimentally at HERA. 




FIG. 3. A graph for 7 + p — » p + jets + X for which diffractive factorization fails because of 
soft color interactions among the spectator partons. 



VI. DIFFRACTIVE PARTON DISTRIBUTIONS FOR (3 -> 

The diffractive parton distribution d A*j x tpi t) I dxjp dt is essentially a long distance 

object, which is not amenable to calculation using perturbative methods. However, Regge 
phenomenology provides an expectation for the behavior of d A*; xjp,t) jdxjp dt when 

the parameter 

(3 = i/xjp (29) 
is small compared to 1. In the case of the ordinary parton distributions, this expectation is 

U/a(x,h) -const, x (l/xf (0) (30) 
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where 5(0) is the pomeron intercept. That is, one expected 

F 2 (x, Q 2 ) ~ const, x (l/xf^ 1 (31) 
for small x. The simplest version of this expectation is that 5(0) is the same as in soft 



pomeron physics, a(0) — 1 ~ 0.08 |15 |. However, this expectation is hard to reconcile with 
evolution: if it holds at some rather small value of Q 2 , then F 2 should be steeper at larger 



values of Q 2 [[RJ. A steeper dependence was also expected on the basis of the perturbative 



version of the pomeron analyzed at leading log level |]r7[. Indeed, a steeper dependence, 



with a(0) — 1 ~ 0.4, is found in experiments at HERA |18| 

The analogous expectation [[| for the f3 dependence of diffractive parton distributions at 
fixed xjp is 

dff^xjp^-x^t) fiy® 

- ~ const, x - . (32) 

dxjpdt \(3) v ' 

In this case, the constant is proportional to the triple pomeron coupling used in Regge 
physics. Presumably, 5(0) here should not be the true pomeron intercept ~ 1.08 but should 
be the larger value, ~ 1.4, found in the HERA experiments. Then, of course, it is debatable 
whether the triple pomeron coupling used here should be the same as that found in soft 
Regge physics. Uncertainty over the precise values, however, should not obscure the simple 
prediction that the diffractive F 2 should behave at small (3 like 

W^EA^ mm ,.JlX m -\ ,33, 
dxjpdt \(3) v ' 

with <j 5(0) — 1 < 0.5. The small (3 behavior of is analyzed in Refs. TTDJ . 



VII. GLUON DISTRIBUTION FOR (3 -» 1 

The diffractive parton distribution d /„/*![(£, /x; xjp, £)/ dxjp dt, like the ordinary parton dis- 
tribution, is essentially not calculable using perturbative methods. Recall, however, that it 
is possible to derive "constituent counting rules" that give predictions for ordinary parton 
distributions f a /A(x,[i) in the limit x — > 1 for not too large values of the scale parameter 



/i [p0|| . In the same spirit, we consider in this section and the next the diffractive parton 
distributions in the limit (3 — > 1 where (3 = £,/xjp. Our analysis is similar to that if Ref. [[3]. 
The authors of that paper concluded, with certain caveats, that the diffractive gluon distri- 
bution should behave like (1 — (3) 1 as (3 — > 1. Our reanalysis suggests a behavior between 
(1 — (3) 1 and (1 — (3)°, depending on how certain nonperturbative issues are resolved. For 



the diffractive quark distribution (not treated in Ref. 0), our analysis in Sec. |VIII| suggests 
a behavior between (1 — (3) 2 and (1 — (3) 1 . 

This analysis involves both hard and soft subprocesses, put together in a manner that is 
not under solid theoretical control. Thus there is not a clearly correct final answer, as far as 
we can see. What we do here is to provide some calculational results that can help restrict 
the range of answers and give some basis for the reader's informed judgment. 
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Define r = (3— 1 = £/xjp — 1. Then we examine d — t)xjp, /i; xjp,t) jdxjp dt in the 

limit r — > 0, after first taking xjp << 1 so as to separate out pomeron exchange from other 
Regge pole exchanges. We consider the diffractive gluon distribution first (a = g). Later, 
we supply the modifications needed for the diffractive quark distribution. Throughout this 
section, except as specifically noted, we consider Feynman graphs in null plane gauge A + = 0. 

Our analysis is based on a model that consists of a selected set of Feynman graphs, 
with the internal loop momenta integrated over a selected integration region. Within this 
model there is a hard subgraph in which all internal propagators are far off-shell and a 
soft subgraph in which the propagators are near to being on-shell. We evaluate the hard 
subgraph at lowest order in perturbation theory. We offer a perturbatively based conjecture 
about the behavior of the soft subgraph in the Regge limit xjp — > 0. 

The simple graphs considered here do not include the graphs that contribute to QCD 
evolution of the parton distributions. Thus we imagine that the analysis applies to the 
diffractive parton distribution at a starting scale /i that is not too large (say, 2 GeV). 
Standard parton evolution starting at this scale will soften the distributions. 



A. Decomposition into hard and soft subgraphs 



According to Eq. (|20|), the diffractive gluon distribution is the square of the matrix 
element of an operator that destroys a gluon with longitudinal momentum fraction (1— t)xjp, 
where the matrix element is taken between the initial proton state and a final state that 
includes the scattered hadron plus anything else. Since a color octet gluon is destroyed, 
while the initial and final hadrons are color singlets, the final state must include at least one 
gluon. We consider here the minimal model in which the final state includes precisely one 
gluon. Then this gluon carries a very small momentum fraction txjp. Call the momentum 
of the final state gluon q^, as depicted in Fig. £|. Since this gluon is on-shell, we have 

* = (34) 

Before proceeding, we pause for a technical point. In general in this section we use the 
null plane gauge A + = 0, but for the final state gluon this is not convenient. The polarization 
vectors e M (g,j) for transverse polarization in the j direction have minus components that 
grow like 1/r in the limit r — > 0: 

c-(?,i) = -T = - ? W- ( 35 ) 
q + XjptPX 

We avoid this singular behavior by changing the polarization vectors for this gluon to A~ = 
gauge. The difference Ae M between the old polarization vector and the new is proportional 
to q 11 , so that changing polarization vectors has no effect after we sum over a gauge invariant 
set of graphs. In A~ = gauge, e~(q,j) = 0, e l (q,j) = 5 % i and 

e + (g,j) = - = ^P/4 (36) 
q~ 
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Thus e M (g, j) is predominantly transverse, with a plus component that vanishes as r — > 0. 

Some vocabulary will be helpful for discussing the physics of time and momentum scales 
in this problem. There are three relevant longitudinal momentum scales. We call partons 
with p + ~ P£ fast partons. For instance, the valence quarks in a hadron are typically fast 
partons. We call partons with p + ~ xjpP^ slow partons. Finally, we call partons with 
p + ~ rxjpP^ very slow partons. The final state gluon is such a parton. 

The gluon cloud surrounding any hadron contains gluons at any momentum fraction. 
Consider a gluon with momentum fraction x and transverse momentum of order m, where 
m « 0.3 GeV gives the scale of a typical hadronic mass or transverse momentum. The 
contribution of such a gluon to the null plane energy P~ of an intermediate state is p~ = 
p 2 /2p + , which is the minus momentum of a free gluon with transverse momentum p and 
plus momentum p + . This "kinetic" minus momentum is of order m 2 /(xP^). Thus a gluon 
of the type that we are calling a slow gluon has a kinetic null plane energy p~ ~ m 2 / {xipP^) 
and survives for a typical null plane time Az + ~ xP^ /m 2 . 

Notice that the final state gluon has a large minus momentum, q~ = q L 2 /(2rxipP A ^), at 
least as long as its transverse momentum q is not too small. Now q is not observed; we 
are to square the matrix element and integrate over q. We cannot say anything about the 
region of very small q, but we can analyze the contribution to the integral from the region 
71 defined by 

q 2 > rm 2 . (37) 

We will consider the contribution to df dlS / dxjp dt from the region 7t with the hope that the 
contribution from the complementary region q 2 < rm 2 , is not large enough to overwhelm 
— or, worse, to cancel — the contribution from region TZ. 

We will also evaluate the contribution from the smaller integration region in which the 
transverse momentum q is large: q 2 ^> m 2 . Since this region is described by rather standard 
short distance dynamics, its contribution should provide a lower bound on the true result. 
(Again, we assume that this contribution is not canceled by some long distance contribution.) 

For q e TZ the minus momentum q~ = q 2 /(2txjpPa) of the final state gluon is large 
compared to the kinetic minus momentum p~ ~ m 2 / Pj[ of a typical fast parton and is also 
large compared to the kinetic minus momentum p~ ~ m 2 / (xipPj[) of a typical slow parton. 
This large minus momentum flows through the graph and is carried out of the graph by 
the detected gluon. Thus the parton detection with r — > creates a hard process that 
happens on a null plane time scale Ax + that is short compared to the typical time scale 
for interactions within the proton or its cloud of slow gluons. We base our analysis on this 
observation, using low order perturbation theory for the nearly instantaneous interaction 
that probes the parton distribution. 

Three graphs that involve the lowest order hard interaction are shown in Fig. [5]. In 
these graphs, the large minus momentum q~ flows through only one or two propagators, as 
indicated by the heavy lines. These propagators are far off shell. We refer to this part of 
the graph as the hard subgraph. The rest of the graph is the soft subgraph. 

There are also two graphs of the same order that involve connecting the g M and gluons 
to a quark line, as depicted in Fig. || These graphs require a separate discussion, which we 
omit in this paper on the grounds that low momentum gluons are likely to be dominant in 
the pomeron compared to low momentum quarks. 
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FIG. 4. Structure of amplitude contributing to the diffractive gluon distribution. At least one 
gluon must be emitted into the final state. 




FIG. 5. Graphs with an order g 2 hard subgraph. 



In Fig. we denote by and s M , respectively, the momenta of the gluon leaving the 
soft subgraph and entering it after the hard interaction. Let the momentum fraction of the 
first gluon be r + /P^ = (1 + a)xjp. Then momentum conservation fixes s + / P + = axjp. 

We integrate over a. It is convenient to distinguish between the two identical gluons by 
requiring that r + > — s + . That is, a > —1/2. In principle, the a integration runs to a = oo, 
but the reader can verify that the region o ^> 1 is not important. We take it as a working 
assumption, to be verified later, that the region \o\ <C 1 is also not important. Thus the 
important integration region is o ~ 1. That is, the gluons that couple the hard subgraph to 
the soft subgraph are typical "slow" gluons. 

We integrate over the transverse momentum r, setting s = r + 1. We also integrate over 
r~, setting s~ = r~ - Mf/(2P/) + (Mf + (P A /) 2 )/(2(l - Xjp)P£). We suppose that the 




FIG. 6. More graphs with an order g 2 hard subgraph. These graphs involve quark exchange 
from the soft subgraph and are not considered in this paper. 
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hadron wave functions fix r and r~ to be no larger than the ordinary size for slow gluons, 
r 2 ~ m 2 and r~ < m 2 /{xjpP^). (Contributions in which r M is part of a hard virtual loop 
are more properly considered to be part of a higher order correction to the hard subgraph.) 

Taking the limit r — > 0, we find that if we consider the hard subgraph to be a function 
of r M , s^, and g M , then it is independent of r~, s~ , r and s and is also also independent of 
Xjp. In addition, the dominant polarizations for the gluons entering the hard subgraph from 
the soft subgraph are transverse; in a shorthand notation, 

J2 [HardHSoftW ~ E [Hard] ^ [Soft] ^ (38) 

H,v={+-,1,2) ij={l,2} 

That the hard subgraph is independent of r~ and s~ is not surprising, since r~ <C q~ and 
s~ < g". That it is independent of xjp follows simply from its invariance under boosts in 
the longitudinal direction. It is not, however, obvious that the hard subgraph is independent 
of r and s and that only transverse polarizations are important. These results follow from 
an argument that we describe in Appendix The picture that emerges is one in which the 
soft gluon cloud surrounding the hadron is probed by an interaction that is effectively local 
in null plane time, x + and in transverse position x. 



B. Structure of the diffractive gluon distribution 

We take advantage of these results by setting r™ = s~ = and r = s = in the 
hard subgraph and restricting the polarization sum to transverse polarizations. Then the 
diffractive distribution function depends on the integral over r~ and r of the soft subgraph. 
We write this structure as depicted in Fig. [7], 



I | roo rod 

^ 2 E Yl j 12 da ' A *ki{v'i x JP^ t; s A , s A >) J i 2 daA i:j (a, x P) t; s A , s A >) 



dxjpdt 64vr 3 r 2 s ^ a , ~ " MV " ' ^' ' 7-i/ 2 

x V f d\ MT n (^cr',rr M% mn (q,a,T) 

Here the first line represents the soft subgraphs while the second line represents the integral 
of the hard subgraphs. We discuss these factors below. 

The function Aij(a, Xp, t; s A , s A r) is the amplitude for the proton to emit a gluon with 
transverse polarization % and momentum fraction (1 + a)xjp then absorb a gluon with trans- 
verse polarization j and momentum fraction axjp. 

A i '(a,xj P ,t;s A ,s A/ )=x IP P+Yl I dy-e^^ v ~ 

x (P A ,,s A ,\T {A{(0, y~, Or) A*(0)} \P A , s A ). (40) 

There is a summation over the color a of the gluon operators; in Ai we average over these 
colors. The factor XjpP^~ arises from changing the integration in Eq. (59) from / dr + to / da. 

We now turn to the hard interaction. We begin by considering the probed gluon, which 
carries momentum According to Eqs. (ET) and (p0|), the operator that probes the gluon 
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distribution is F^ m (y). (Only v e {1,2} contributes in Eq. (f20|).) In A + = gauge, this is 
simply d + A™(y). (We have changed the polarization vector e(q,n) for the final state gluon 
from e + = 0, so this gluon could couple to the A + A m operator in F^ m (y). However, this 
graph is not allowed when the two exchanged gluons are in a color singlet state.) The d + 
becomes a k + = (1 — t)xjpP^~ ~ xjpPa that we absorb into the normalization. This leaves 
a propagator for the probed gluon, 

—iD mil (k) —i —i 

" {) =T^-J2^f)e^k,f) = -^-e»(k,m). (41) 
k z + ie k z + %ty k z + xe 

Thus each probe operator gives a 1/k 2 factor and leaves the amputated hard interaction 
graph dotted into a polarization vector for a gluon of momentum polarization m, in 
A + = gauge. 

We call the amputated hard interaction graph ^^""(q, er, r). The indices m and n are 
the transverse polarizations of the probed and final state gluons respectively; a and h are 
their colors. Recall that the final state gluon polarization is represented by a polarization 
vector e tl (q,n) is, by our convention, in A~ = gauge, as in Eq. fl3~6[) . In A4, the exchanged 
gluons are approximated as having zero transverse and minus momenta; one integrates over 
these momenta in the definition of A. The exchanged gluons have transverse polarizations 
i and j respectively. 

In the r — > limit, the factor k 2 in the denominator is 

k 2 ~ -q 2 /r. (42) 

Thus the probed gluon is far off shell, k 2 3> m 2 , as long as q is in the integration region 
Tl, q 2 > Tin . Similarly, the virtual lines internal to A4 are far off shell and have a simple 
form when q G 71. This is so even though the region TZ includes small transverse momenta, 
q 2 m 2 . Indeed, the small transverse momenta near the boundary of TZ dominate the 
integral in Eq. fl39D , as we will see. 

Next, we examine the soft subgraphs and then the hard subgraphs in some detail. 



C. The soft subgraph 

In this subsection, we investigate the xjp dependence of A^(a, t; sa, sa>) in the limit 
xjp — > 0. Comparing the expected Regge form (|8|) of df^/dxip dt with Eq. (0), we see that 
the expected Regge form of A^ is 

A 13 (a, xjp, t; s A , s A >) oc %>° W . (43) 

The question is, what is the pomeron trajectory a(t) within the context of the analysis used 
in this paper? 

The operators Eq. (|40|) are in A + = gauge, but the amplitude can be put in a gauge 
invariant form by reexpressing it in terms of the gluon field operators F^ u defined in Eq. (|18|) . 
We simply replace 

s+A> . d+A* . F+i 
A 3 -> — -> -i— -> -i— 44 

s + + te s + + te s + + te 
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R (l+cr')x. 



FIG. 7. Factored structure of the diffractive gluon distribution function for r — ► 1. 

and make a similar replacement for A . The ie choice here is of some significance. We will 
discuss it in the following subsection. The resulting form for A^ is 



A lJ (a,xjp,t;s A ,s AI 



1 



Efdy 



-^icrxpP+y 



xjp(g + ie)(l + a)P A H „ 

x <P A ,, s A ,|T {F+ J (0, y', T )F+\0)} \P A , s A ). 



(45) 



Since the propagators in the soft subgraph are not far off-shell, one is not really justified to 
use perturbation theory to investigate A^(a, xjp, t\ s A , s A r). Nevertheless, perturbation theory 
is suggestive. Consider Eq. (|^) using Feynman gauge. Suppose that the gluon annihilated 
by the operator F+ l connects to an on-shell quark carrying no transverse momentum and 
plus momentum p + = XP A , with A of order 1. That is, the gluon couples to a "fast" quark. 
The relevant factor is 



IE 



U(\'P£,s)igt c ^ 



A ■ 



S). 



(46) 



where A' = A — £jp(l + a). In the limit xjp — > 0, this is independent of xjp. Similarly, the 
coupling of the operator to a fast quark gives no xjp dependence. Thus perturbation 
theory suggests that the operator matrix element in Eq. fl4"5] ) is independent of xjp for small 
xjp. Considering that there is a factor 1/xjp in Eq. (f45D , we find that the pomeron trajectory 
appearing in Eq. (f43D is a(t) = 1 in this most naive analysis of the soft subgraph. This is 
close to the pomeron trajectory observed in nature, but of course this most naive analysis is 
too naive, and we expect that soft interactions among the gluons modify the result. What 
we see here is that the picture of the pomeron as two gluon exchange, which often gives 
results that are surprisingly good considering the simplicity of the picture , works rather 
well also in this context. 
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D. Small a singularity 



The analysis of the preceding subsection has revealed a 1/a singularity within the in- 
tegration domain — l/2<o~<ooof the momentum fraction a. What is the nature of 
this singularity? If we stick to A + = gauge, one can check that it arises from the 
singularity in the gluon propagator, 



s 2 + it 



(47) 



This singularity is usually interpreted with a principle value prescription, but there is no 
compelling reason for this choice. In fact, A + = gauge is not an effective tool for examining 
the nature of the 1/a singularity, since in this gauge the singularity is a gauge artifact. Thus 
we choose to examine this question in Feynman gauge, in the style of Ref. fl~2f . (Arguably, 
it would be better to do the whole problem in Feynman gauge, but this leads to its own 
complications.) 




FIG. 8. Graph with a singularity near s + = 0. 

We consider the graph shown in Fig. |] in Feynman gauge but with a transverse polar- 
ization chosen for the gluon carrying momentum r M . It is helpful to choose a frame in which 
xjpPj^ ~ m. Then 



P+ = (1 - XIP )P+ > m P A ,= Pl,/[2(1 - xjp)P+] < m 
q + = TXjpP^ <tim q~ = /[2rxjpP^] > m. 



The graph contains the structure 



s 2 + it (q+ s) 2 + it 
The denominator (q + s) 2 has the form 



(48) 



(49) 
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(q + s) 2 + ie = 2q- s + s 2 + ie = 2(q~ + s~)s + - q- s + 2g + s - s 2 + ie (50) 

Now q~ is large and positive while s~ is small in the dominant integration region. Thus the 
denominator has the approximate form 

1 



2q-[s+ + --- + ie] (51) 

where the dots indicate small terms. This denominator contains the only important depen- 
dence on s + as long as \a\ < 1 in s + = axjpP^ . Notice that there is a singularity very 
near to s + = 0, but that, since there are no other singularities nearby, we can deform the 
integration contour away from this singularity. Thus we deform the s + contour into the 
upper half complex s + plane. On the deformed contour, we have s + ~ m or o ~ 1. 

Now we examine the numerators. The largest component of the quark current J M is 
J + , while the largest component of the current N u of the final state gluon is N~ . Thus 
J M iV M ~ J + N~ . Thus the dominant term in J^iV 74 is obtained by replacing 

W -> ^ s^N". (52) 

The factor s^N^ is approximately s + N~ as long as s + is not small, and we know that s + 
is not small since we have deformed the integration contour so that it does not approach 
s+ = 0. 

The next step is to restore the s + integration contour to the real axis, taking care not 
to cross any singularities. This means that we should move the l/s + singularity in Eq. ([52] ) 
infinitesimally into the lower half s + plane, so that our replacement becomes 

JpN* -> — T — s^N". (53) 

This is in keeping with the usual notation in which integration contours are along the real 
axis, with poles infinitesimally displaced from the integration contour. 

The replacement (|53|) gives the dominant contribution to our graph. However, if we 
attach the gluon carrying momentum s M everywhere in the hard subgraph and sum the 
leading terms obtained with this replacement, we will get zero because of the Ward identities 
obeyed by the hard graph and because the two gluons carrying momenta and s M together 
form a color singlet. (The relevant identities are discussed in the appendix of Ref. jT2f ). 

We have thus encountered the bete noir of Feynman gauge: the leading contributions 
graph by graph come from unphysical polarizations that cancel when one sums over graphs. 
What we need is the subleading contribution. That is easy. We replace 



J»N» = — J — s,N» + J,N». (54) 
s + + ie 



J" = ; . 55 

s+ + ie v ' 



where 



Now we throw away the first term in (Q), since it will cancel, and keep the remainder. 
This gives the structure in Eq. (|44[), but now with a prescription for the l/s + singularity 
determined on physical grounds. 
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E. The hard subgraph 



We now turn to the hard interaction function M.. It is a simple matter to evaluate this 
function. We find 



4(a + ie)(l + a 



5, 



lib 



a(l + a) 6 ij 8 mn - a 5 in 5 jm + (1 + a) S im S jn 



(56) 



Thus M. oc r° as r — > 0. The ie prescription in the 1/a factor arises from Eq. fl5~ll) . It matches 
the ie prescription in the soft function Aij, Eq. ([15]) , so that, after a contour deformation, a 
is never much smaller than 1. 

Inserting Eq. ( pUf) and A; 2 ~ — q 2 /r into Eq. (p^|), we obtain 



c?Xip (it 647r 3 r 2 _ 



x 



/DO roo 
da'A* kl (a',xip,t- : s A ,s A >) / daA ij (a,x IP ,t; s A , s A >) H khj '(a, a') 
-1/2 J-1/2 

d\ fl(q 2 > tM 2 ^ 



X i(2vr)2 (cf/r) 2 " (57) 

Here Ti 3 (a, a') is a rational function of a and a' that is not of particular interest. Recall 
that we integrate q over the region TZ defined by q 2 ^> rm 2 (Eq. fl3T|)), where m is a typical 
hadronic mass or transverse momentum. Here we make the prescription more precise by 
integrating over q 2 > tM 2 where M is any fixed mass such that M ^> m. 



F. Result 



Performing the integration in Eq. (|57|) gives for the r — > limit of the diffractive gluon 
distribution 

d fg/A _ 1 1 ^ ^ 



dxjp dt 256tt 4 M 2 2 , 



X 



/OO /"OO 
dcr'A^^XiPjtjSAjSA') / daA ij {a,x IP ,t;s A ,SA')'H kh:i ((y,(7'). (58) 
-1/2 J-1/2 



Note that dfgj A /dxipdt is independent of r as r — > 0: 

dxjpdt \ hj h ^ v > 

This may seem surprising, since the distribution of gluons in a typical hadron behaves like 
fg/A^x; fio) ~ const, x (1 — x) p with a rather high power p « 5. 

The leading behavior comes from the lower endpoint of the integration in Eq. fl57|). Thus 
it is sensitive to the cutoff chosen. Recall that we took q 2 > tM 2 because as long as 
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q 2 /r ^> m 2 , the internal lines in the hard subdiagram are far off-shell. This appears to 
be the natural cutoff. As soon as the internal lines of the subdiagram through which g M 
flows are not far off-shell, the the gluon that goes into the final state and the probed gluon 
can attach at different space-time points in the gluon cloud of the hadron. Then there is 
the opportunity for cancellation, as both gluons sample the color charge of the hadron as 
a whole and find that the hadron as a whole is a color singlet. It is difficult to check this 
conjecture directly in a realistic model of nonperturbative hadron structure. However we 
have checked in a very simple model where all the graphs can be included exactly. This 
model is too simple to have the correct pomeron behavior, but we find that it does have r° 
behavior in the r —>■ limit. The model is described in Appendix [B]. 

One might reasonably conjecture that a larger cutoff would be imposed by nonper- 
turbative physics in a more realistic model. For instance, if the gluon emitted into the 
final state effectively had a substantial mass m g , then q~ would be (q 2 + m 2 ) / '(2xjprP A + ). 
This would induce an effective cutoff q 2 > m 2 in Eq. (p7|). Then we would have obtained 
dfgj^(Pxjp, /i; xjp, t)/dxjp dt oc (1 — f3) p with p = 1. This corresponds to the style of analysis 
of the constituent counting rules and gives the result (1 — f3) 1 found in Ref. ||. Presumably, 
the contribution from q 2 3> m? must be present and is not likely to be canceled, so that 
df^(Pxjp, /i; xjp, t)/dxjp dt should not be smaller than (1 — (3) 1 at large (3. That is, the power 
p should not be larger than 1. 

We conclude that if the diffractive gluon distribution is parameterized as 



dfg/AifalP'KXlP't) 



oc(l-/3) p (60) 



dxjp dt 

for (3 — > 1 at moderate values of the scale /i, say 2 GeV, then 

< p < 1. (61) 

The choice p ~ corresponds to an effectively massless final state gluon, while p ~ 1 
corresponds to an effective gluon mass. 



VIII. QUARK DISTRIBUTION FOR (3 -> 1 

In analogy with the gluon case, we write 



dxjpdt 64tt 3 t 2 7-i/2 J (2tt) 2 J -1/2 J (2vr 
xS* a (a', s'; xjp, t; s A , s A/ )S^(a, s; xjp, t; s A , s A >) 
d 2 q 



x EE 

SqS k I J 



(27T) 



2 



MutSL k , t] s q , s k )* MZish k > s , °, T ', s q , s k ) , . 

(k 2 - it) (k 2 + it) ' 1 ' 

Here the soft function S^a, s; xjp, t; s A , s A >) is the amplitude for the proton to emit a gluon 
with polarization fi, momentum fraction (1 + a)xip and transverse momentum r = s — 
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Pa', then absorb a gluon with transverse polarization v and momentum fraction axjp and 
transverse momentum s: 

S^{a, s; xjp, t; s A , s A >) = xjpP+ £ / " dy~ / dy e"*^^^-^ 

x (P A ,, ^'|^(0, y" y) ^(0)|P A , s A ). (63) 



The function S 1 is simply related to the operator matrix element A l \ Eq. fl40|), that appeared 
in our discussion of the diffractive gluon distribution: 

ds 

S ll {a, s; xjp, t; s A , s A >) = A tJ (a, xjp, t; s A , s A >). (64) 



(2k) 



The function 



M^(cL^8,(r,T;s q ,8 k ) (65) 



represents the amputated hard interaction graph. Here q is the momentum of the antiquark 
that enters the final state and k is the transverse momentum of the probed quark. We have 
k = — q — 1 by momentum conservation. The variables and s q are the helicities of the 
probed and final state quarks respectively; I and J are their colors. 

There is an important difference with the gluon case. The quark has a mass m q . Thus 
the minus momentum of the on-shell quark entering the final state is 

Then 

y? « (q2 2 + r '"' ) (67) 

in the r — > limit. What counts here is the mass of the final state quark as it emerges from 
the hard interaction and propagates into the final state. Presumably the best model for m q 
in this role is the constituent quark mass (~ 0.3 GeV), not the much smaller current quark 
mass. This is a substantial mass, so that the condition that defined whether the virtual lines 
in M. are far off shell, 



(q 2 + ml) 



T 



> m\ (68) 



is satisfied for any q 2 when r is small. Thus we do not need to restrict the integration 
region. On the other hand, the region q 2 <C m 2 q is not important in the integration. 

We evaluate M. in the r — > limit using null plane spin defined with the plus direction 
as special for the probed quark and defined with the minus direction as special for the 
final-state antiquark. We find 



M"" x^fsPSu f = w(s^T^w(-s q ). (69) 
/ q^ + m l 
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The w(s) are two component spinors w(+±) = (1, 0) and w(— §) = (0, 1). Then we can write 
the 2x2 matrix T using transverse Pauli spin matrices a 1 and a 2 . For transverse indices 
[iv, we find 

1 . , 1 



a~(s ■ o)<j 3 + a 3 (r ■ a)a l 



a + it 1 + a 

+2<F(k • a) + — — crY — a j q l 

a + it 1 + o" 

+if — <rV + oV + 2^ m q (70) 

V a + ze 1 + a J q 

(We hope that the Pauli spin matrices a 1 will not be confused with the momentum fraction 
a that occurs in this equation in the combinations 1/a and 1/(1 + a).) For one transverse 
index and one plus index we find 

T +3 « 2^' (71) 

Also T* + ph r +? , while T ++ does not give leading contributions as r — > 0. Finally, we note 
that r _1/ and are not needed because they multiply in A + = gauge. 

The spin function T^ u is rather complicated, but we are concerned with only two of its 
properties. First, if we think of V in coordinate space as a function of the separation 
between the points where the two exchanged gluons attach, then T is proportional to S(y + ) 
and to a linear combination of 5(|/t) and d5{yx) / dy^- Thus the interaction is hard in the 
sense of being local in y + and yx- Second, and most important, V is independent of r. 

We can now insert Eqs. ( |67D and ( j69|) into Eqs. (|62"D to obtain the r dependence of the 
diffractive quark distribution: 

<m _ 2 J_3Cll_ 1 r r d 2 s r d 2 s> 

dxjpdt T 64vr 3 16 2 s £f Ai J a ° J (2tt) 2 J (2vr) 2 

xS* a (a', s'; Xip, t; s A} SA')S^(a, s; xn>, t; s A} s A ') 

f d\ Tr{r^K,s^,q)tr^( ( T,s,q)} 



(2tt) 2 (q 2 + m. 



2^)3 



2 



The crucial feature here is the factor of r 

We conclude that the constituent counting result for the diffractive distribution of quarks 

is 

*^loc(l-ff. (73) 
dxjp dt 

However, suppose that we interpret the calculation of the previous section as saying that 
the diffractive distribution of gluons is proportional to (1 — (3)° for (3 near 1 when the scale 
/i is not too large. Then the evolution equation for the diffractive parton distributions will 
give a quark distribution that behaves like 

^^" .(l-H) 1 . (74) 
dxjp dt 
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when the scale \i is large enough that some gluon to quark evolution has occurred, but not 
so large that effective power p in (1 — f3) p for the gluon distribution has evolved substantially 
from p = 0. A signature of this phenomenon is that the diffractive quark distribution will 
be growing as \x increases at large (3, rather than shrinking. Perhaps this is seen in the data 



IX. CONCLUSION 

We close with some observations concerning the implications for HERA physics of the 
discussion presented here. 

We have discussed two kinds of factorization relevant to diffractive hard scattering. Both 
are experimentally testable. We say that diffractive factorization holds if the cross section 
is a standard partonic hard scattering cross section convoluted with a diffractive parton 
distribution. The diffractive parton distribution gives the distribution of partons in the 
hadron under the condition that the hadron is diffractively scattered. (If there is a second 
hadron in the initial state and we do not demand that it be diffractively scattered, then the 
cross section should also contain a convolution with the ordinary parton distribution within 
that hadron.) We say that Regge factorization holds if the diffractive parton distribution 
is a product of xfp? a ^ times a pomeron-hadron coupling times a function f a /p((3, t, /f) that 
is interpreted as the distribution of partons "in" the pomeron. These properties together 
constitute the Ingelman-Schlein model 0. 

Diffractive factorization is something that can be disproved for a given process by a 
counterexample at some fixed order of perturbation theory (with wave functions for the 
hadronic states). Correspondingly, it could in principle be proved to hold at any fixed order 
of perturbation theory, which one would take as a strong indication that it holds beyond 
perturbation theory for that process. Diffractive factorization makes no statement about 
the Regge structure of the nonperturbative factors. Regge factorization does make such a 
statement. It would be interesting to study Regge factorization from the point of view of the 
BFKL pomeron, perhaps extracting a model for the distribution of partons in the pomeron. 

Let us discuss first diffractive deeply inelastic scattering, which is simpler than diffractive 
hard scattering processes with two hadrons in the initial state. For diffractive deeply inelastic 
scattering, we argue that diffractive factorization is a consequence of perturbative QCD, 
although a detailed proof is beyond the scope of this paper. From the point of view of 
current theory, Regge factorization for the diffractive parton distributions is a conjecture 
based on experience with soft diffractive scattering. 

The HERA experiments |IJ have now provided evidence concerning these issues. It is a 
consequence of diffractive factorization that the diffractive structure function should 
exhibit approximate scaling as Q 2 is increased with fixed xjp and (3. This is confirmed by the 
data. The dependence on xjp in the form x]p 2a ^ predicted by Regge factorization, Eq. (|i2|), 
is also confirmed by the data. 

For the future, it will be important to measure the diffractive parton distributions in as 
complete detail as possible, using charged and neutral current events, F^ lfT and F^ lS , and 



probes for heavy flavors in the final state p2 |. Especially crucial is the diffractive gluon 



distribution. This can be measured using diffractive deeply inelastic scattering with high 
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Pt jets detected in the final state. If we demand that we see two jets instead of the usual 
single struck quark jet, and if these two jets have a high transverse momentum relative to 
the direction determined by the sum of their momenta, then the hard process is of order 
ota s instead of just a. For such a process, gluons participate as initial partons on the same 
footing as quarks. Since the diffractive quark distributions are already known, it should be 
possible to extract the diffractive gluon distributions. 

Our analysis suggests that the diffractive gluon distribution is quite hard, with a behavior 
between (1 — /3) 1 and (1 — /3) for large j3 = ^jxjp at moderate values of the scaling parameter 
fi, say 2 GeV. The corresponding behavior of the diffractive quark distribution, which is quite 
directly measured in F^, is between (1 — (3) 2 and (1 — (3) 1 . Here the (1 — (3) 1 for quarks 
would arise if the diffractive gluon distribution is large and behaves like (1 — /3)°, so that 
the quarks at large (3 are produced by g — > q + q. 

Given a complete set of diffractive parton distributions, it will be interesting to test the 
evolution equation (|2lf). 

Let us now turn to hard processes with two hadrons in the initial state, as exemplified 
by 7 + P —> jets + p + X at HERA, where we look at the hadronic part of the photon. In 
this case, the factorization that holds in inclusive hard scattering is expected to break down 
in diffractive hard scattering, as shown by counterexamples at a fixed order of perturbation 
theory fH],|5| . If one extracts diffractive parton distribution functions from deeply inelastic 
scattering and uses them to predict cross sections for 7 + p — > jets + p + X, then the 
observed cross section should contain extra terms that do not match the prediction ||. In 
particular, there should be extra contributions that correspond to the jets carrying almost all 
of the longitudinal momentum of the pomeron. Perhaps this corresponds to the "superhard" 
component seen in the UA8 experiment 0] in p + p p + jets + X. 

In summary, the HERA results JIJ, together with the earlier UA8 results jjJJ, have con- 
firmed the basic features of the Ingelman-Schlein picture of diffractive deeply inelastic scat- 
tering. The experiments have shown that diffractive scattering is related to exchanges of 
quanta that, when examined with a hard probe, appear to be the pointlike quarks and 
gluons of QCD. Much more remains to be done, but already we are challenged to connect 
the theory of pointlike gluons to the soft color dynamics that is presumably responsible for 
diffractive scattering. 
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APPENDIX A: STRUCTURE OF THE HARD SUBGRAPH 

In this appendix we investigate the structure of the hard subgraph for the diffractive 
distribution of gluons in hadron in the limit £/xzp — > 1. Recall that the hard subgraph 
is a function of momenta r M , s M , and g M , where r M and s M are the momenta of the gluons 
exchanged with the soft subgraph and q^ is the momentum of the gluon that goes into the 
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final state. Since q^q^ = we will consider the hard subgraph to be a function of q and q~, 
replacing q + by q 2 / (2g _ ). The momentum of the detected gluon is = r M — s M — g M . 
In Sec. |V111| , we studied the hard subgraph in the limit that applies when (,/xjp — > 1. 



This limit is really a dual limit. We take r 



xjpPt but 



+ 



2s + q~ 



< 1. 



(Al) 



We also recall the definition ([37]) of the integration region considered for q, 



q ^ 2 

— ^> m . 

T 



(A2) 



That is, q ^> m? / (2xjpPj[) . We combine this with the assumption that, in the effective 
integration region for r, s, r~ and s~, these variables have values typical for slow gluons: 



m 



r ~ s ~ 



2 2 2 

r ~ s ~ m . 



(A3) 



Then 



r 

— «1 ; 
9~ 



2s+q~ 



«1, 



2s+g~ 



< 1. 



(A4) 



We claimed in Sec. |V11 A| that in this limit 
variables r, s, r~, and s 



the hard subgraph is independent of the 



s, r~, and s~ . We also claimed that the transverse components of the hard 
subgraph dominate over other components in the limit considered, as in Eq. (^). In this 
appendix, we substantiate these claims. 

We consider first the question of the independence on the variables r, s, r~, and s~, 
taking, for the moment, only the transverse components of the hard scattering subgraph. 
This is the same as multiplying the hard subgraph by purely transverse polarization vectors 
for the two gluons exchanged with the soft subgraph. 



We consider the hard amplitude A4, defined as in Sec. [VII E[ Thus Ai does not include 
the propagator for the detected gluon with momentum but does include an e + = gauge 
polarization vector for this gluon. It also includes an e~ = gauge polarization vector for 
the gluon with momentum g M that enters the final state. 

As a matter of convenience, we will analyze the first graph in Fig. |^. Essentially the same 
argument covers the second graph, while the third graph, with a four gluon interaction, is 
quite trivial. 

The Feynman rules give Al as a rational function of the components of r M , s M , and q^. 
We are interested in particular in the behavior of A4 for small values of r~ , s~ , r, s and q 
as specified in Eqs. (|A1|) and (|A4j). For this reason, we need to know if there are any factors 
of the small variables in the denominator of A4. 

The gluon propagator in Fig. || is 



iD^(q + s) 
(q + sf 



+ q + ) + (s + qY5 u _ + 5^(s + q) 1 



{s+ + q+) [2(s+ + q+)(q~ + s~) - (s + q) 



(A5) 
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Setting q + = q 2 /(2g ), the denominator becomes 



2( S +)V 1 + 



q 2 \ J" / q 2 \ / s \ (s + q) 2 



2s + q~ I IV 2s+g J \ q~ J 2s + q~ 



(A6) 



Then using Eqs. ( |AI| ) and (|A4[ ) (in either order), the denominator becomes 

2( S +)V (A7) 

That is, the denominator does not contain factors of the small variables, so that it has a 
finite limit as the small variables tend to zero. 
We can write Ai as a product 

M ijmn = e( S ,z)ae(r,j)^e(A:,m) 7 e( g ,r i ) 5 M^' 5 . (A8) 

The transverse components of the polarization vectors have the form 

e(p,i) I = 5 lI (A9) 

for p 11 stands for any of the momenta r M , s^, g M or k^. For p 11 = r M or we are (tem- 
porarily) defining e(p, i) + = e(p, i)~ = 0. The polarization vector for the detected gluon has 
e(k,m) + = but has a non-zero minus component 

e{k,m)- = ?- = r + \ \ . (A10) 
k + r + — s + — q + 

The polarization vector for the final-state gluon, which we take to be in e~ = gauge 
according to Eq. fl36]) , has a non-zero plus component 

n 

e(g,n) + = ^-. (All) 
Q 

Thus neither the gluon propagator nor any of the polarization vectors contains a factor of a 
small variable in the denominator. 

We now consider the transverse tensor J\A. % i mn as a function of the variables 

{r + ,s + ,r",s",g",r,s,q}. (A12) 

M. is a rational function of these arguments and has dimension D = and boost dimension 
B — 0, where B gives the scaling under boosts in the z direction. That is, for any four- vector 

p + has D = l, B = 1, 
|p| has D = 1, B = 0, 

p- has D = 1, B = -1. (A13) 

Since has B = and D = 0, it can be written as a function of a reduced number of 
arguments, each of which has B = D = 0: 
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We are interested in the limit in which the last five arguments are small, and we know that 
none of these arguments occur as factors in the denominators. Thus when Eqs. and 
hold, M. approaches the limiting value 

M ijmn { T —i 0, 0,0,0,0 j . (A15) 



We thus confirm the claim made in Sec. |VII A| that we can neglect r, s, r , and s in 



M.. Furthermore, since the limiting function is covariant under rotations about the z axis, it 
must be a linear combination of the tensors 8^8 mn , 5 m & n and <5 m 5- ?m . Each of these tensors 
multiplies a rational function of r + /s + = (1 + er)/cr. That is, the coefficient functions are 
rational functions of a. This is, of course, just the structure given in Eq. fl56l) . The argument 
given above does not establish that the limiting function is nonzero, but this is what we 
found by calculation. 

We now return to the issue of whether the transverse components of the hard subgraph 
dominate over other components, as claimed in Sec. |VII A| . We consider 



M af3mn = e(k, m) 7 e(q, n) s M a ^ s (A16) 

and choose values other than {1,2} for a or (3 or both. Now a = — or (3 = — are not 
possible: M a P mn multiplies the soft subgraph with corresponding indices, call it S a/ 3, which 
vanishes for a = — or (3 = — because of the gauge condition A + (x) = (that is, A-(x) = 0). 
Thus we should consider a = + or (3 = +. Let us consider a = + with (3 = j G {1,2} 
as an example that illustrates the general argument. Thus we wish to investigate whether 
M +jmn S +j is dominated by A^' mn <% in the limit specified by Eqs. O and (T£% We 
need an order of magnitude estimate for the ratio of S + j to <Sy . An analysis of lowest order 
graphs indicates that this is the same as the ratio of the corresponding components of the 
polarization vectors for a typical slow gluon, e~ (s^ 1 , i) / (s^ , i) ~ m/ (xjpP^). We write this 
as 



Sij 



(A17) 



and compare 



s 



s 



M +jmn (A18) 



to j\A. % i mn . The analysis is simple. J\/[ + i mn has dimension D = and boost dimension 



B = 1. Thus it can be written as \Js + /q times a function J\f^ mn that has dimension 
D = and boost dimension B = and, like J\A^ mn , has no factors of the small variables in 
its denominator. As our previous analysis shows, M^ mn has a finite limit as the the small 
variables tend to zero. (Actually, J\[^ mn vanishes in this limit because of its spin structure, 
but we will not need to use this fact.) We note that the factor that multiplies J\fi mn 

(A19) 




y/s + q~ 

vanishes in the limit specified by Eq. (|A4|) . This establishes the claim. 
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APPENDIX B: A MODEL 



In this appendix we compute the diffractive gluon distribution df^/dxjp dt in a simple 

model. As we will see, this model does not exhibit the pomeron behavior x^ a ^ with 
a(t) « 1. Nevertheless, the model provides a check that df^ A (Pxjp;xjp,t)/dxjpdt can have 
a finite limit as (3 — > 1 in an exact numerical evaluation of the graphs to a given order of 
perturbation theory. 

The model we will use is scalar-quark QCD, given by the Lagrangian, 



2.-7.. X na ni» L f^\2 1 fa A ^2, f Faddeev-Popov ] 

1 



C = D^qD^q - m 2 qq G a nv G^ g A (qq) (d ■ A a Y + < } 



+-(^) 2 --M 2 2 -G#g 

= ^QCD + £<j>- (Bl) 

This model has soft and collinear singularities just as in QCD. As one further simplification, 
we model the diffractively scattered particle by the scalar <ft field with a (f>qq interaction to 
the quarks. Then the perturbative <pqq interaction plays the role of the nonperturbative 
Bethe-Salpeter wave function of a real QCD meson. The probability for finding a gg-pair 
inside the meson in this model falls of as 1/k 4 in the ultraviolet, where k is the transverse 
momentum of the quark. This good behavior in the ultraviolet is due to the fact that G car- 
ries the dimension of mass. Thus the model allows for a simple treatment that approximately 
simulates properties of bound quarks inside a hadron. 

Using this model, we compute the diffractive gluon distribution function within a meson, 
dfg/A.(P x iPi xp,t)/dx]pdt, as given in Eq. (|22"D , to the lowest nontrivial order in g, order g A . 
One of the order g A diagrams for the function G^ A (P, q, k) in Eq. ( p2|) is shown in Fig. One 
gluon is detected, and a single gluon is exchanged between the two quark loops on opposite 
sides of the final state cut. This gluon carries that portion of the momentum transfer from 
the meson that is lost to the detected gluon (and in a physical process would be lost to 
the hard interaction). The heavy bar ending with gluon lines represents the gluon density 
operator in Eq. (|T^). In addition to the graph shown, at each loop the gluons can also 
attach to the lower quark line. Also in both cases there is a graph where the gluon lines are 
crossed. In addition one gluon can attach on each of the quark lines. Finally there are the 
two-gluon two-quark contact interaction graphs. In total this implies 8 2 = 64 combinations. 
After symmetry considerations, there are four types of amplitudes that must be explicitly 
evaluated. The four amplitudes are shown in Fig. |TU|. 

The final form for G^ A that we obtain is, 

G ^ (P ' g ' k) ~ 2(2^)3/5(1-/3)^ J d q (FTiM (B2) 

where k + = j3xipP + , 

c ^ k ) = T^Tf ( fc+ V + ^ W - W + 9t - Kk + g+), (B3) 
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j 

FIG. 9. Diagram contributing to the diffractive gluon distribution in a meson in the model of 
this appendix. 




and G Mi/ (P, q, k) is a second rank tensor containing the sum of quark loops with all possible 
combinations of attachments by two gluons. 

We have evaluated the loop integrals in Eq. (|B^) in terms of their explicit Spence function 
expressions. For this we have used the computer code of G. J. van Oldenborgh and J. A. 



M. Vermaseren [p3 |. Their method is an independently formulated extension of the Form 
algorithm of G. 't Hooft and M. Veltman [|24[| . The authors of |23| claim that their algorithms 
provide easier isolation of both asymptotic behavior and potential numerical instabilities. 
In one of the most sensitive regions, that of low t, their algorithms have proven to be more 
accurate. 

To summarize the specifications of our calculation, it was for small t, one external line 
was massless, and one had to integrate over a large region in the transverse space of the 
variable q. Due to gauge invariance, several constraints are placed on which we have 
numerically checked to hold. G tiV {P, q, k) has the form 



CT(P 1 q t k) = y £ l Tf'Aj 
J=l 



(B4) 



Here the Tj are tensors made from the available vectors P, q, k. The Aj are scalar functions 
of the dot products of these momenta. In general there would be ten terms in the sum, but 
gauge invariance cuts this in half. To find the possible Tj, define 



q 2 P^-P-qq", 
k 2 P v -P-kk\ 
k- qP^ - P ■ qk^ 
S u k =q-kP u -P-kq u . 



DM 

g 

pv 

CM 



(B5) 



Notice that q ■ P q = 0, etc. Now define 



rpHV _ pfj, pi/ 
1 l — r q r k 1 
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FIG. 10. Four graphs for the amplitude in Fig. 
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J 3 


— D q Jr k 5 


rpflU 


— o ? o fc , 


rpflV 

J 5 


= q-kgT 



(B6) 



These have the property that they are polynomials in the components of the momenta and 
that q^Tj V = and k u Tj V = 0. Thus G^ v is properly gauge invariant. In our calculation, 
we transformed G^ v into a basis where the first five tensor forms are those given above and 
the other five are 



rpflU 
rpflU 

J 10 



r k r q i 
pM 

r k °q 5 
1=1 k r q ' 
D k D q 5 



(B7) 



We have checked over a wide range of the parameters xp, t and the transverse momentum 
q that the coefficients of the tensors in Eq. ( |B7| ) vanish, leaving only the tensors in Eq. ( |B6| ) 
with nonvanishing coefficients. As a second check for gauge invariance of G^, we have 
numerically checked that q^G^ u = and k v G^ v = 0. 

In Fig. |ll| we show the diffractive gluon distribution function multiplied by fixp, at 
\t\ = 1 GeV 2 , and for a selected choice of xjp values. The masses of the quarks and mesons 
were M 2 = m 2 = 0.1 GeV 2 , while the couplings were G = 0.3 GeV and g = 0.1. We see 



that the model does not exhibit the behavior 
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df*f A ((3xjp;xjp,t) 2 

— ^— j -r. oc Xjp (B8) 

axjp at 

for xjp — > 1 at fixed (3 that would be characteristic of pomeron exchange. One needs at least 
one rung of a gluon ladder to obtain this behavior. The model exhibits the behavior 

J g/A VH IP IP, J _ o 

as /3 — > 1 at fixed xjp, as in Eq. (|59|). 




FIG. 11. Diffractive gluon distribution 0xjpdfg^(^X]p;x]p,t)/dxipdt in the model of this 
appendix as a function of /3 for (from top to bottom) ayp=0.01, 0.02, 0.05, 0.1 with \t\ = 1 GeV 
In this model, there is no dependence on a renormalization scale /i. The parameter choices are 
M 2 = m 2 = 0.1 GeV 2 , G = 0.3 GeV and g = 0.1 

One can understand the (3 and xjp behavior seen in this numerical study from an analytic 
viewpoint. In Sec. [VII| , the detected and emitted gluons coupled do a slow gluon. In the 



present simple model, they couple directly to a fast quark. This means that the virtual 
quark line in Fig. |9] is far off shell when q 2 /((l — /3)xzp) 3> m 2 . For example, it is off shell 
for q 2 ~ m 2 for any (3 as long as we take xjp 1. This gives quite a different behavior from 
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that found in Sec. |VII] . In the model, the loop serves as a built-in low momentum cut-off. 
However, the meson can be replaced by a single fast quark line if we substitute an artificial 
low momentum cut-off 

q 2 > (1 - /3)xip m 2 (BIO) 
Then a simple power counting analysis gives, 

dfif A ^Xjp- xjp, t) const. [°° , 2 9(q 2 > (1 - (3) Xjp m 2 ) 

d( lT r 2 TTT-, m iTo 



dxjp dt (I - (5)x 2 P Jo 1 {ql/[{l - 0) xp}} 

for xjp 1 and (1 — j3) 1. Performing the integral gives 



dxp dt 



(Bi2) 



This agrees with our numerical findings. We find it reassuring that a fully consistent field 
theoretic calculation of the diffractive gluon distribution gives results that agree with expec- 
tations similar to our analytic arguments in Sec. [VI]. 
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